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Using a Luttinger-liquid approach we study the quantum fluctuations of a Bose-Josephson junc- 
tion, consisting of a Bose gas confined to a quasi one-dimensional ring trap which contains a localized 
repulsive potential barrier. For an infinite barrier we study the one-particle and two-particle static 
correlation functions. For the one-body density-matrix we obtain different power-law decays depend- 
ing on the location of the probe points with respect to the position of the barrier. This quasi-long 
range order can be experimentally probed in principle using an interference measurement. The 
corresponding momentum distribution at small momenta is also shown to be afi'ected by the pres- 
ence of the barrier and to display the universal power-law behavior expected for an interacting ID 
fiuid. We also evaluate the particle density profile, and by comparing with the exact results in the 
Tonks-Girardeau limit we fix the nonuniversal parameters of the Luttinger-liquid theory. Once the 
parameters are determined from one-body properties, we evaluate the density-density correlation 
function, finding a remarkable agreement between the Luttinger liquid predictions and the exact 
result in the Tonks-Girardeau limit, even at the length scale of the Friedel-like oscillations which 
characterize the behavior of the density-density correlation function at intermediate distance. Fi- 
nally, for a large but finite barrier we use the one-body correlation function to estimate the effect 
of quantum fluctuations on the renormalization of the barrier height, finding a reduction of the 
effective Josephson coupling energy, which depends on the length of the ring and on the interaction 
strength. 
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I. INTRODUCTION 

The possibility to study Bose-Einstein condensates 
confined to ring traps constitutes one of the frontiers of 
the experimental progress with ultracold atomic gases [l| . 
The nontrivial topology of these traps together with the 
uniformity of the potential along the ring circumference 
makes them an ideal system for investigating persistent 
currents and superfluid properties of the gas. 

While current experimental setups display relatively 
weak transverse (i.e. radial) confinements, we analyze in 
this work the case where the transverse confinement is 
so strong that only longitudinal (i.e. tangential) quasi- 
ID motion is allowed along the ring. In a quasi-lD ge- 
ometry, the phase coherence properties of the gas are 
drastically changed with respect to their 3D counter- 
parts. Phase fluctuations destroy true long-range or- 
der, and by increasing the interaction strength the gas 
changes from a quasi-condensate, i.e. a condensate with 
fluctuating phase [3, Q to a Tonks-Girardeau gas 
where repulsions are so strong that they mimic the ef- 
fect of Pauli pressure in a Fermi gas and the conden- 
sate is strongly depleted. Here we consider such a one- 
dimensional ring trap containing a localized repulsive po- 
tential which creates a "weak link" connecting the two 
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ends of the loop (see Fig. [T]), a situation that may be 
viewed as a realization of a Bose-Josephson junction. 
Bose-Josephson junctions have been already experimen- 
tally realized using a double- well geometry and arrays [5| . 
In the configuration considered here, quantum fluctua- 
tions tend to destroy the phase coherence along the ring, 
while the tunneling of bosons between the ends of the 
loop favors a well-defined phase difference across the bar- 
rier. We will study the interplay between these compet- 
ing effects. 

We start by investigating how the presence of the bar- 
rier affects the quantum fluctuations and hence the coher- 
ence properties of a Bose-Josephson junction in a quasi- 
lD ring at arbitrary values of the interaction strength. 
In the absence of the barrier, these properties have been 
extensively studied, employing a variety of techniques, 
from low-energy Luttinger-liquid approaches and confor- 
mal field theory [aHK to exact methods in the Tonks- 
Girardeau regime [3, Q ■ In the presence of an infinitely 
high barrier, using the Luttinger liquid approach for a 
finite ring, we evaluate here the first-order correlation 
function, which describes the decay of phase coherence 
along the ring, and recover a previous result from con- 
formal field theory Q. We then compare the results 
for the one-body density matrix and for the particle- 
density profile along the ring with the corresponding 
exact results in the Tonks-Girardeau limit of impene- 
trable bosons obtained through a Bose-Fermi mapping 
method This enables us to determine the numerical 
values of the nonuniversal parameters of the Luttinger- 
liquid theory. Knowing these parameters we estimate 
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FIG. 1: Sketch of the Bose-Josephson junction on a ring trap 
studied in this work. 



two-particle properties such as the density-density corre- 
lation function. Finally, we turn to the case of a large but 
finite barrier, treating the tunneling across the barrier as 
a perturbation. We use the results for the first-order cor- 
relation function in the infinite-barrier limit to estimate 
the effect of the quantum fluctuations on the effective 
height of the barrier, i.e. on the Josephson coupling en- 
ergy. As a result, we predict how the renormalization of 
the Josephson energy depends on the ring length and on 
the interaction strength. 



II. LUTTINGER LIQUID DESCRIPTION FOR A 
RING WITH A LOCALIZED BARRIER 



We consider N bosons confined to a uniform, quasi-lD 
ring-shaped trap of circumference L to which a local- 
ized repulsive potential Vbarr(a;), located at a; = = L, 
has been superimposed. The bosons interact with each 
other through a repulsive contact potential v{x — x') = 
g5{x — x'). The corresponding Hamiltonian in terms of 
the bosonic field operators ^'(a;), ^'^(a;) reads 



n 



2m 



dx vl/t (x) (a;) * (a;) * (a;) . 



(1) 



The presence on the ring of the (very large) barrier po- 
tential will be taken into account below by imposing open 
boundary conditions at a; = = L (see Sec. Ill Bp and by 
adding a tunnel term (see Sec. IVip . 
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The parameters K and Vg are related to the microscopic 
interaction parameter of the original Hamiltonian ([T|) , 
the phase field (j){x) is related to the velocity of the fluid 
v{x) = hV(t){x)/m and the field 9{x) defines the fluctua- 
tions in the density profile p{x) according to 



p{x) = [po + n(a:)] 
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where po = N/L is the average density of the fluid, 
n(a:) = \'9{x)/Tr and 9b is a constant fixing the posi- 
tion of the first particle with respect to the origin of the 
X axis. The fields 11 and 4> satisfy the commutation rela- 
tion [i, [H 



[n{x),<j){x')] ^iS{x-x'). 



(4) 



Note that this approach includes not only the lowest- 
order hydrodynamic expression for the density fluctua- 
tion at long wavelength: the higher order terms in the 
sum ([U enable the description of the discrete nature of 
the particles up to distances a 1/po- Our approach 
does not allow however to probe shorter distance scales 
because of the assumption of linear phonon dispersion 
modes in the effective Hamiltonian ([5]). The bosonic field 
operator is obtained from the hydrodynamic expression 
^-+(3;) v/p(a:)e-''^(^) and reads 



'f^x) = A^ypo + n(x) 



, (5) 



where ^ is a nonuniversal constant, the value of which 
depends on the way the Luttinger liquid approach is reg- 
ularized in the short wavelength limit. This issue will 
be discussed in Sec. IV Al below, where the value of the 
constants A and 9b will be fixed. The field operators ^' 
and as defined through Eq. ([5]) satisfy the standard 
bosonic commutation relations as a consequence of the 
commutation relation ^ among the field operators H 
and d. 



A. Low-energy theory 

In order to evaluate the equilibrium correlation func- 
tions at large and intermediate distances we adopt the 
Luttinger liquid approach, i.e. we approximate the sys- 
tem Hamiltonian ([1]) by the following effective low- 
energy Hamiltonian in terms of the fields 9{x) and 0(a;) 
which describe the density and phase fluctuations on the 



B. Mode expansion of the Luttinger fields 9 and ((> 
with open boundary conditions 

We start by considering the case of an infinitely high 
barrier, which corresponds to a ring with open boundary 
conditions. In Sec. IVII we will treat the case of a large, 
finite barrier by considering the tunneling among the two 
sides of the barrier as a perturbation. 
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In order to evaluate the first- and second-order corre- 
lation functions for the bosons on the ring junction, we 
derive here the expansion of the fields Il{x) = V9{x)/tt 
and (f>{x) in terms of canonical bosonic annihilation and 
creation operators bk and &^ satisfying the commutation 



relations 6^,] 



Specifically, we expand the op- 



erators (t){x) and n(x) in Fourier modes for x £ [0, L] 



-t-oo 



</>(x) = ['^l,.K■+'^2,,<]e''=^■^ (6) 

+ 00 

U{x) = J2 [ni,,6fc,+n2,,6ije''=^-, (7) 



j = -aD 



where we have set kj — ^pj- The constant p and the 
complex coefficients '/>ij, 4'2.j: ^i.j find Tl2.j are deter- 
mined by imposing three constraints: (i) the commu- 
tations rules [n(x),0(x')] = iS{x — x'); (ii) the open 
boundary conditions, which imply vanishing current den- 
sity i.e. V0(O) = V0(L) = 0, Vn(0) = Vn(i) = 0; 
(iii) reduction of the Hamiltonian to the diagonal form 



In order to take into account 

that we are using an approximate, long wavelength the- 
ory, we introduce a short distance cutoff a ^ pq~^ in the 
sum over the modes. The final result reads 



(x) = (t>o 
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— ^ cos{Trjx/L) e 
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n{x) = n, 







^/j cos{tt jx/L) e 21, 



bk,+b 



the latter implying 



(x) ~ ttIIoX 



,(8) 



(9) 



—p sm{njx/L)e 2t 
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III. EXACT DESCRIPTION IN THE 
TONKS-GIRARDEAU LIMIT K = 1 

In the limit of infinitely strong repulsion between the 
bosons, which corresponds to the value if = 1 for the 
Luttinger liquid parameter, an exact solution exists for 
the bosonic many-body wavefunction $(a;i, ...xat) (in 
first quantization). We shall use it throughout this pa- 
per in order to test the results of the Luttinger-liquid 
theory in the limit K — \^ thereby fixing the values of its 
nonuniversal parameters. 

The solution, due to Girardeau is obtained by 
mapping the bosons onto a gas of noninteracting, spin- 
polarized fermions subject to the same external poten- 
tial. The bosonic many-body wavefunction ...xat) 
is then obtained in terms of the fcrmionic one as 

...Xjq) = ...X]^)^f{xx...X]^), (11) 

where the mapping function ^(xi, ...xat) — 
ni<j<^<7vsign(a;j —xi) ensures the proper symmetry un- 
der exchange of two bosons. The fermionic wavefunction 
is given by $_f(xi, ...xat) = (l/VM) det[-0j(xfc)]j,fc=i...Ar, 
ipj{x) being the single particle orbitals for the given 
external potential. Note that <& p vanishes every time two 
particles meet as required by Pauli's principle, and hence 
describes well the impenetrability condition g oo for 
the bosons. In our specific case, the orbitals for a ring of 
circumference L and open boundary conditions are 



iPj{x) = y/i2/L) sm{TTjx/L) 



(12) 



with j = 1, CXI. 

As a consequence of the Bose-Fermi mapping, all the 
bosonic properties which do not depend on the sign of the 
many-body wavefunction coincide with the correspond- 
ing ones of the mapped Fermi gas. This is the case 
e.g. for the particle density profile and for the density- 
density correlation function. Other properties like the 
one-body density matrix and the momentum distribution 
are instead markedly different for bosons as compared to 
fermions. In particular, the calculation of the one-body 
density matrix requires in principle the calculation of a 
(N-l)-dimensional integral, which is known to simplify in 
some cases. Examples are the homogeneous gas with pe- 
riodic boundary conditions Q or the case of a harmonic 
confinement [9]. 



where kj = nj/L {p = 1/2). The zero mode IIq is directly 
related to the particle number operator through normal- 
ization: Ho = [N - No)/L where Nq = {N) = poL. It 
is conjugate to the zero-mode phase operator (/>o such 
that [no,(/'o] = i- Using this fact, one can explicitly 
check the commutation rule between 9(x) and (l){x') from 
the mode expansions ^ and pH)) : it turns out to be 
[9{x), ^(x')] = i7ru(a; — x'), where u is the unit step func- 
tion, consistent with (HJ. Finally, inspection of the di- 
agonalized form of the Hamiltonian yields the linear dis- 
persion relation Ukj = hvskj for the modes. 



IV. ONE-BODY DENSITY MATRIX AND 
MOMENTUM DISTRIBUTION IN THE 
INFINITE-BARRIER LIMIT 

In this section we focus on the one-body density ma- 
trix and on the momentum distribution for the case of a 
bosonic ring of circumference L with an infinite barrier, 
with the aim of analyzing the differences with respect to 
the case of an infinite system, as well as to the case of a 
ring in the absence of the barrier. 
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A. Contribution from phase fluctuations to the 
one-body density matrix 

The one-body density matrix, defined as G{x,x') = 
(\I>^(a;)5'(a;')) yields a measure of the coherence along 
the ring. It is possible to measure the one-body density 
matrix and off-diagonal long range order experimentally 
by measuring the interference pattern of atomic matter 
waves coming from two holes in the trap (see e.g. 'll] for 
the case of a cigar-shaped 3D Bose gas). According to 
Eq. ([5]) the bosonic field operator has three contributions: 
(i) the phase (j){x) (ii) the density fluctuation n(a;), and 
(iii) the higher harmonics of order 2m9{x) of the density. 
The most important contribution to the one-body density 
matrix at large distances is the one due to the phase fluc- 
tuations which correspond to the lowest-energy modes of 
the bosonic fluid in the ring (see e.g. 0, 01)j while the 
two latter contributions give rise to subleading correc- 
tions which we do not analyze further here. 

To lowest order we approximate the bosonic field 



operator ([5]) as ^{x) 



A^fp^e 



the prob- 



lem then reduces to the computation of the quan- 



-i4>{x) i4>{x')\ 



Since 



tum average Qq{x,x') — A^Pq{ 
the Luttinger-liquid Hamiltonian Q is quadratic in 
the field 0(a;) we immediately obtain Qq{x,x') — 

^^pocxp^— i (^[4i{x) — (/)(a;')]^^^ . The phase-phase cor- 
relation function is evaluated with the help of the 
mode expansion ([8]); using the fact that the ground- 
state average over the bosonic modes is ((fej. + &[)(&; + 
bj)} = Ski and the property T.'jLi } ^'"^ coa{jj) = 
— i In [l — 2 cos 7 e~" -I- e~^"] one readily obtains 

{cp{x)cj^{x'))^-^\n[{n/Lf {a' + d\x - x'\2L)) 

X {a^ +(f{x + x'\2L))\ , (13) 

d{x\L) — L\sv[\{ttx/L)\/tt being the cord function. This 
leads to 



Qo{x,x') = pobo^Q 

P^VK + d^{'2x\2L)] [q2 + d^{2x'\2L)] 
+ d'^{x - x'\2L)] [a2 + d'^{x + x'\2L)] 



(14) 



where we have introduced the nonuniversal constant 
^0,0 = \A.\^{poa)^ . The above expression (fT4|) yields 
the leading-order term for the one-body density matrix 
at large distances. By taking the limit a — > we recover 
the result obtained in [6(| using the methods of conformal 
field theory. 

If the distance among x and x' is large compared to the 
cutoff length a, the one-body density matrix displays a 
power-law decay of the form Qq(x,x') oc \x — x'\~'~', where 
the exponent 7 can be derived from the expression (|14p . 
and in particular depends on the location of the probed 
points [yj. Indeed, if the two points are away from the 



obtained in the thermodynamic limit [l3| whereas if they 



approach the edges {i.e. x a and L 



r' < 



a) the 



exponent is 7 = a result known in the context of 
quantum phase fluctuations in a ID superconducting wire 
of length L 12] . In the case where one point is at one edge 
and the other in the bulk we obtain 7 — . These three 
different behaviors are illustrated in Fig. [2 where we plot 
the one-body density matrix Goix, a;') as a function of x' 
for various choices of the probe point x. In the same 
figure we display also the behavior for a homogeneous 
ring in absence of the barrier, obtained by a procedure 
analogous to the one outlined above. 



C(x,x')=PoVo 



a'^ + d'^{x-x'\L) 



(15) 



Note that, as the coordinate x' runs along the ring, in the 
presence of the barrier the coherence decreases monoton- 
ically, while if the barrier is absent coherence is recovered 
as x' approaches L — x. 

The different power-law behaviors are in principle ob- 
servable for a quasi-lD Bose gas in a ring trap geometry; 
it is required to have a high barrier well localized on a 
length scale a. 



B. Momentum distribution 

We proceed by studying the momentum distribution 
n{q), obtained by Fourier transformation of the one-body 
density matrix with respect to the relative variable, 



n{q) = 



dx 



da;'e-*«(^-^')g(x,x'). (16) 



edges one finds 7 



2if ' 



which corresponds to the result 



We have resorted to a numerical calculation for the eval- 
uation of the momentum distribution taking as input the 
one-body density matrix obtained in Eq. p^ . This al- 
lows to estimate the main features of the momentum 
distribution at wavevectors q smaller than the cutoff 
wavevector qc 1/a. The behavior at large wavevectors 
q ^ qc needs an accurate treatment of the short-distance 
behavior of the many-body wavefunction [l^, [l3| and 
is beyond the regime of validity of the Luttinger-liquid 
method. The result for the momentum distribution is 
illustrated in Fig. [3] for two values of the boson num- 
ber in the ring, and at varying interaction strength. As 
a general feature (see the inset of Fig. [3]), we observe 
that at intermediate values of q the momentum distri- 
bution displays a power-law behavior n{q) ~ q'^/i^^)-^ 
with the same power predicted for a homogeneous ring 
(see e.g. [3])- This result is readily understood as the 
different power laws described in Sec. lIVAl only occur at 
the edge of the integration region with a negligible weight 
with respect to the bulk contribution. Still, by compar- 
ing the details of the momentum distribution of the ring 
with the barrier with the momentum distribution of a 
uniform ring, (see the main panel of Fig. [3]), we find that 
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FIG. 2; (Color online) Top and middle panels: one-body den- 
sity matrix in logarithmic scale (arbitrary units) as a function 
of the coordinate x' (in units of the ring circumference L) for 
various choices of the coordinate x with respect to the posi- 
tion of the barrier, located at a; = 0: top panel a; ~ 0, middle 
panel x = L/2. Bottom panel: one-body density matrix in 
logarithmic scale (arbitrary units) taken at a:: ~ and a;' ~ L 
as a function of the length L of the ring (in units of the average 
interparticle distance a). In each panel we plot three values 
of the parameter K (from bottom to top K = 1, K = 2 and 
K = A); the solid lines correspond to the results from Eq. (|14|) 
and the dashed lines are the solution (llSp for a homogeneous 
ring (periodic boundary conditions). The dotted line in the 
top panel is the exact solution for a Tonks-Girardeau gas in 
the thermodynamic limit. The linear behavior of t/o in log- 
arithmic scale correspond to the predicted power law decays 
with various exponents 7 (top panel, 7 = 1/2K, middle panel 
7 = i/AK and bottom panel 7 = 1/K). 



in the presence of the barrier the momentum distribution 
is decreased at small momenta. This is in agreement with 
the fact that the barrier reduces the coherence along the 
ring. The result is reminiscent of the one obtained for a 
ID gas in presence of disorder [l^, where the reduction 
of the momentum distribution at small momenta is also 
observed. 



q/qc 1 

FIG. 3: (Color online) Momentum distribution n{q) in units 
of 27r|^pQ as a function of the wavevector q in units of the 
cutoff momentum qc — l/a for = 10 bosons on a ring with 
an infinite barrier (solid lines) and for a homogeneous ring in 
absence of the barrier (dashed lines) for various values of the 
Luttinger parameter K (from top to bottom K = 4, 2, 1). 
The inset shows the same quantity (in logarithmic scale, ar- 
bitrary units, same line conventions as the main figure) eval- 
uated for = 10^ bosons. The dotted lines indicate the 
predicted power law decays q^^^^'^^~^. 



V. PARTICLE-DENSITY PROFILE AND 
DENSITY-DENSITY CORRELATION 
FUNCTION IN THE INFINITE-BARRIER LIMIT 

Extending the quantum average techniques outlined in 
Sec. llVl to the limit of an infinitely high barrier it is pos- 
sible to evaluate also the inhomogeneous particle-density 
profile and the density-density correlation functions. In- 
terference between particles incident on and reflected by 
the barrier leads to the occurrence of Friedel-like oscil- 
lations in the density profile and in its correlator, which 
are typical of strongly correlated ID fluids. We describe 
here these oscillations within the Luttinger liquid ap- 
proach, for any value of the coupling strength, finding 
that they are more and more marked as the coupling 
strength increases. In the Tonks-Girardeau limit of in- 
finite boson-boson repulsion we compare the predictions 
of the Luttinger-liquid approach with the exact results, 
which enables us to fix the nonuniversal parameters of 
the latter. 



A. Friedel oscillations in the particle density profile 

We compute the particle-density profile by taking the 
quantum average {p{x)) of the density operator ^ on 
the ground state, namely 

m— — 00 

(17) 

To evaluate the quantum averages we exploit the fact 
that the Hamiltonian ^ is quadratic in the field 9{x) 
(and recall that n(a;) — V9{x)/ti). We use the property 



6 



that if X and Z are Gaussian variables then (^X e^ 

(XZ)e3<^'>. Hence, we are left with the evaluation of 
the {9{x)9{x')) correlation function, to be taken at equal 
points X — x'\ the (Yi{x)9{x')) correlation function is ob- 
tained from the previous one by taking the derivative 
with respect to the variable x. Using the mode expan- 
sion (llOp and a procedure similar to the one outlined in 
Sec. II VI we obtain 



{e{x)e{x')) = 



'k'^{UI)xx' 
q — in 



x'\2L) 



a"^ + d'^{x + x'\2L) 



(18) 



and 



{n{x)e{x')) = n(nl)x' + 

K 



K 



d{x + x'\L) 



27r a"^ + <P{x + x'\2L) 
d{x - x'\L) 



2tt a^ + d^{x - x'\2L) 



sign(x — a;')(19) 



The final expression for the density profile reads 



(P(^)> 



Po 



= 1 + 2E 



-'^ \ d{2x\2L) 



mK 
t^Pq 



sm(2rmT pqx + 2m9B) 



cos{2rmTpox + 2to0b) 
d{2x\L) 



d'^{2x\2L) 



(20) 



The density profile is modulated by oscillations with 
wavevector multiples of 2npQ. Notice that in the case 
K — 1^ where the system can be mapped onto a non- 
interacting spin-polarized Fermi gas, the wavevectors of 
the oscillation are multiples of 2kp^ where kp = npo is 
the Fermi wavevector, and hence correspond to the well- 
known Friedel oscillations For the case of generic 
K the m = 1 oscillations decay with the power law x~^ 
(see e.g. 

Let us now concentrate on the case K — 1. In the 
thermodynamic limit (L — > cxd, N oo, at fixed po = 
N/L) the expression (PU)) for the density profile at short 
distances (to 0{l/x)) reduces to 



{P{^)) 

PO 



1 



a cos{2'Kpf)X + 26b) 



(21) 



This can be compared with the thermodynamic limit of 
the exact expression derived using the Bose-Fermi map- 
ping I, [3 



PO 



1 - 



sin(27rpo2;) 
2t:pox 



(22) 



allowing us to fix the coefficients a and 6b to the val- 
ues a — l/(27rpo) and 9b = ir/i. Note that the lat- 
ter choice for 6b is in agreement with the condition 
7^ 0,±7r,±27r. . . obtained by imposing that the par- 
ticle density profile should vanish at cc = and x = L Q . 
Once the constants a and 9b are chosen, the constant A 
in Eq. ([5]) can be fixed by comparing the expression for 




FIG. 4: (Color online) Particle density profile {p{x)) in units 
of the average particle density po as a function of the spatial 
coordinate x along the ring (in units of the ring circumference 
L) for various values of the parameter K. Main figure, K = 1 
(solid line: result from the Luttinger-liquid model, dotted line, 
exact result from the Bose-Fermi mapping); inset, K = 2 (top 
green line) and K — 4 (bottom blue line). 



the coefficient 6oo entering Eq. for the one-body den- 
sity matrix with the exact value 6q™'^* = 2~^/'^y/7reAQ^ ^ 
0.521 f§, ^ where Ag = 1.282 ... is Glaisher's constant. 
The resuh is \A\^ = 21/6,^61/2^^6 ^ ^ rpj^-g ^^^^^g 

has been used in plotting Fig. [21 

In Figure |3] we illustrate the density profiles for vari- 
ous values of K, obtained by the Luttinger-liquid expres- 
sion ((20)) using the above choice for a and 9b d9']. The 
figure displays also the exact result for the density pro- 
file obtained from the Bose-Fermi mapping, (x) = 
J2i=i where the single-particle orbitals ipj{x) 



are defined in Eq. p2)) . The comparison shows how 
our choice of parameters a and 6b reproduces extremely 
well the density profile oscillations even on a finite ring. 
The figure also illustrates how the Friedel oscillations dis- 
play maximal amplitude in the strongly interacting limit 
K ^1. 



B. Density-density correlation function from 
Luttinger-liquid theory 



We turn now to the density-density correlation func- 
tion 5'(a;, a;') = {p{x)p{x')) — {p{x)) {p{x')) . This quantity 
encodes the information on the structure of the fiuid, 
i.e. on the correlations between density modulations at 
different parts of the fluid, while it vanishes for an ideal 
Bose gas. The Fourier transform of the density-density 
correlation function with respect to the relative variable 
is directly accessible experimentally by light-scattering 
methods (see e.g. [20j and references therein). 

The density-density correlation function is obtained 
with the quantum average method described in Sec. IIVI 
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andlVl One has to compute 

+ OC 



{p{x)p{x')) 



E 



g2i(m— m )dB ^i27rpo{mx — m x ) 



m,m — — oo 



X ([po' + Po(n(a;) + n{x')) + n(a;)n(x')] 



and HH) is 



(n(x)n(x')>=-2^ 



(P{x - x'\2L) - cos(7r(x - x')/L) 



{cP{x~x'\2L) + a'^Y 
(f{x + x'\2L) - cos(7r(x + x')/L) 



(n 



The average can be performed using the gen- 
eral result for Gaussian variables (^XYc^'^ = 

((Xr) + (XZ)(rZ))e^<^'>. The novel correlator 
needed for the calculation in addition to Eqs. pS]) The final result reads 



{d^{x + x'\2L) + a^f 



(24) 



+ 00 



^2i(m — m')dB p*27rpo {mx — m'x') 



{p(x)p{x')) = (Poa)('"'+'""^^e2-- „.,..g.. 

X [p2 ^_ (n(a;)n(x')) + 22po((n(a;) + n(a;'))(m6l(a;) - m'e{x'))) 
~A{Il{x){me{x) - m'6i(a;')))(n(a;')(m6i(a;) - m'e{x'))) 

mm' K / — 9 \ k'm'^K 



+ (f{x + x'\2L) 
a"^ + <P{x- x'\2L) 



a"^ +(P{2x\2L) 
I 



+ (P{2x'\2L) 



Sm'^K 



(25) 



This equation displays the general structure of the 
density-density correlations to all orders in m and m', 
and by considering only the first terms of the expansion 
m, m' = 0, ±1 we recover the known results [1, 0]. 

We proceed by comparing the density-density correla- 
tion function S{x,x') with the exact result for K = 1. 
The latter is obtained from the Bose-Fermi mapping 
as [13: 



5"= 



[X, X ) 



N 



-I 2 



J = l 



(26) 



where the single-particle orbitals ^j{x) are defined in 
Eq. (dll). 



Figure [H] displays the results obtained from the 
Luttinger-liquid method at various values of the Lut- 
tinger parameter using the choice of parameters a 
and 9b determined from the density profile in Sec. IV Al 
and compares to the exact ones in the case K = 1. The 
agreement found is very good, even for the Friedel-like 
oscillations at wave vector k ~ 27rpo; this is at the bound- 
ary of the expected regime of validity of the Luttinger- 
liquid theory and illustrates how a reasonable choice of 
the non-universal parameters in the effective model al- 
lows for surprisingly accurate predictions. 




-0.05 



FIG. 5: (Color online) Density-density correlation function 
S{x,x') (in units of po) from the Luttinger-liquid model as a 
function of the coordinate x' (in units of the ring circumfer- 
ence L), with a; = L/2 and for various values of the Luttinger 
parameter K {K = 1, top red line, K = 2, middle green line, 
K = 4 bottom blue line). The dotted line corresponds to the 
exact solution obtained from the Bose-Fermi mapping in the 
case K — 1. 



VI. RENORMALIZATION OF JOSEPHSON 
ENERGY BY QUANTUM FLUCTUATIONS 

In this section we consider the effect of a finite barrier 
on the ring located at x = = L, and described by a 
localized potential Vbair(a;) = Td{x), where T is the tun- 
neling strength of the barrier. This yields the following 
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barrier term in the Hamiltonian ([T]): 

= r*t(L)^(0) + h.c. (27) 

The above equation takes into account the possibihty for 
a boson to tunnel through the barrier potential, and h.c. 
refers to the Hermitian conjugate corresponding to tun- 
neling events in the opposite direction. In the hydrody- 
namic formulation for '^^x), the transfer of one boson is 
ensured by the operator exp(— it/?) where ip — (j)(L) — (f)(0). 
By neglecting the density fluctuations in the field opera- 
tor ([5]) we then recover the usual Josephson Hamiltonian 

Hj = Ej cos ^, (28) 

where Ej = 2p()T is the Josephson energy of the junc- 
tion. 

Quantum fluctuations of the bosons in the ring on 
both sides of the barrier tend to smear the phase ip 
and hence suppress the tunneling strength. Indeed, from 
the diagonal Hamiltonian of section HlBl the ring consti- 
tutes an oscillator bath for the junction with linear spec- 
trum fkuk = T^hvsj /L; the resulting model is very similar 
to the one describing a superconducting Josephson junc- 
tion coupled to a resistive environment [2l|. Tunneling 
events thus induce excitations of the modes of the ring 
with energy between Hluq = ttHvs / L and the high energy 
cutoff huJii = nhvs/ct. 

When the Josephson energy is smaller than the low- 
est mode hujQ, corresponding to small rings L < L* ^ 
hvs/Ej, the junction can be treated as a perturbation 
and every mode modifies Ej. The effective Josephson 
energy results from averaging Tij, Eq. (^5)) . with respect 
to the ground state of the unperturbed Hamiltonian ([2]): 
Ef = {Uj) with {n.j)/Ej = go(i,a)/po. Then 

Ef = Ej[^) fori<r. (29) 

The Josephson energy decreases with the power 
law L^^/^ of the one-body density for probe points at 
the edges of the ring. This case includes the limit of an 
infinitely high barrier, where Ej —f 0, L* — )■ oo, and 
which is illustrated in Fig. [21 bottom panel. 

When Ej is larger than huiQ (with Ej < hujh), only 
the modes with energies larger than Ej contribute to the 
renormalization. Consequently, to obtain the effective 
Josephson energy we need to average over wavelengths 
between a and the characteristic length £ = ttHvs/Ej^: 
Ef = ((Hj)) with 

Ef/Ej = GoiL, a)/go{e, a) ^ {a/iy/'' (30) 

The effective Josephson energy in this case is obtained 
by solving the above self-consistent equation with respect 
to £, with the result 

Ef^Ej^^^j iorL>L\ (31) 




L/q 200 



FIG. 6: (Color online) Renormalized Josephson energy as a 
function of the length of the ring. From top to bottom K = 
4, 2, 1 with hujh = IQEj. Ef decreases as and 
reaches a constant value at L = L*. 

In this case Ef is independent of the ring circumfer- 
ence L. Our results ([29l) and ((3T|) are summarized in 
Fig. [51 As a main conclusion, we find that quantum fluc- 
tuations dramatically reduce the tunnel amplitude with 
respect to its bare value entering the Hamiltonian (|27p . 
especially in the case K = \ . Note however that the re- 
duction saturates at a nonzero level for rings larger than 
the so-called healing length L* . The continuity between 
the two regimes L < L* and L > L* defines the healing 
length of the ring 

As a final remark, we would like to mention that our 
approach is equivalent to the renormalization group for- 
malism [2^ 1 or the self-consistent harmonic approxima- 
tion [m. 

VII. SUMMARY AND CONCLUDING 
REMARKS 

In summary, in this paper we have studied the equilib- 
rium properties of a quasi-lD interacting Bose gas con- 
fined in a ring trap with a localized barrier. In the limit 
of infinite barrier we have studied the coherence, den- 
sity profiles and density-density correlations of the gas 
using a Luttinger-liquid approach and the quantum av- 
erage method. Our results recover and extend those pre- 
viously known by the use of conformal field theory meth- 
ods. As physical consequences of our analysis, we find 
that the one-body density matrix, when probed at var- 
ious points with respect to the barrier position, is ex- 
pected to display universal power law behaviors with dif- 
ferent exponents which depend only on the Luttinger pa- 
rameter K. We also find that our method permits to 
describe accurately the Friedel oscillations (due to the 
presence of the barrier) occurring in the particle density 
profile and in the density-density correlation function. 
Once the non-universal parameters entering the effective 
model are fixed by comparing the density profile to the 
exact one in the Tonks- Girardeau case K = 1, we find 
that the Luttinger liquid model well agrees with the ex- 
act result for the density-density correlation function at 
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a length scale which is at the boundary of the validity of 
the Luttinger-liquid model. 

The analysis performed in the limit of infinite barrier is 
then use to study perturbatively the presence of a large, 
finite barrier. By taking into account the effect of quan- 
tum fluctuations we find that the effective Josephson en- 
ergy (i.e. the tunnel amplitude across the barrier) is 
reduced with respect to its bare value, in a way which 
depends on the length of the ring with respect to a typ- 
ical healing length, a maximal reduction occurring for 
long rings. 

The effect of the renormalization of the tunnel am- 
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